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The strong interaction effects of isospin- and spin-dependent 
nucleon-nucleon correlations observed in many-body calculations are in- 
terpreted in terms of a one-pion exchange mechanism. Including such 
effects in computations of nuclear parity violating effects leads to en- 
hancements of about 10%. A larger effect arises from the one-boson 
exchange nature of the parity non-conserving nucleon-nucleon interac- 
tion, which depends on both weak and strong meson-nucleon coupling 
constants. Using values of the latter that are constrained by nucleon- 
nucleon phase shifts leads to enhancements of parity violation by factors 
close to two. Thus much of previously noticed discrepancies between 
weak coupling constants extracted from different experiments can be 
removed. 







The problem of determining the parity-violating interaction between nucleons has 
drawn recent attention with measurements of the 133 Cs anapole moment [1], the TRI- 
UMF measurement of proton-proton scattering [2] and the interactions of epithermal 
neutrons with heavy nuclei [3]. The parity non- conserving NN potential, V PNC , of 
Desplanques, Donoghue and Holstein (DDH) [4], constructed so that the only un- 
known quantities are the weak meson-nucleon coupling constants, which govern the 
strength of the weak meson-nucleon vertex function, has been the standard tool to 
analyze these data because the parity-violating observables have been expressed as 
linear combinations of these, presumably fundamental, weak coupling constants. As 
explained in various reviews [5-9], a problem occurs because different values of the 
weak coupling constants are needed to describe different experiments. In particular, 
the 133 Cs anapole moment requires a pion-nucleon weak coupling constant /„- that is 
larger than that predicted in Ref. [4], but the observation of circularly polarized 7 
rays in the decay of 18 F data require a very small value for Furthermore, a recent 
analysis of heavy compound nuclei [10] finds parity violating effects that are larger 
(by factors ~ 1.7 — 3) than those predicted using DDH potential. 

We next explain why the results of nuclear structure calculations [11,12] cause 
us to examine the effects of spin-isospin, nucleon-nucleon, correlations on parity vi- 
olation. Nuclear parity violating effects have been typically analyzed in terms of 
a parity-violating single nucleon shell model potential, constructed from the DDH 
potential using a Hartree-Fock approximation and RPA correlations (e.g. [5,10,13]). 
However, two-particle-two-hole correlations are known to be as vital elements of nu- 
clear structure. The spin-independent effects of the short range correlations are often 
incorporated [5] using the Miller- Spencer correlation function [14], but this does not 
take into account all of the correlation effects [15,7]. 

In particular, recent variational studies of nuclear matter [11] and 16 [12] have 
demonstrated that spin-isospin correlations are very important. To be specific, con- 
sider two-nucleon pair correlation functions defined by the expectation values 

Pr(r) = -L^ttl £tf(r - r tJ )F h ^}, (1) 

i<j 

with = |rj — and F^j are various two-nucleon operators. For the central term 
T c = 1. The work of Refs. [11,12] found that with F it j = t,- l ■ Tj cri ■ crj, (F = or) 
or Fij = Ti ■ TjSij (r = tenr), the pr(r) can be relatively large, even larger than the 
well-known effects of short-distance repulsion. Other operators F cause much smaller 
effects and are ignored here. 

What is the impact of spin-isospin correlations on calculations of parity-violating 
observables? Consider the construction of the single-particle PNC potential fj PNC 
which governs the interaction between a valence nucleon (i) with a spin + , N = Z 
core. In using the Hartree-Fock approximation, one folds the operator \/ PNC with the 
density matrix of the core. The pion exchange term contains the operator i(rj x Tj) z , 
where j represents the core nucleons. The expectation value of this operator vanishes. 
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But if one includes spin-isospin correlations, the relevant matrix elements include a 
factor Ti-Tj, and the product x Tj) z Ti-Tj = —2(Ti — Tj) z —i(Ti x Tj) z contains 
a non- vanishing term, —2t{. Thus a new non- vanishing contribution will appear. 

The starting point for our evaluation of the PNC single-particle potential is the 
effective parity- violating nucleon-nucleon potential [4] between two nucleons (i,j): 



~9 P [h p r t ■ r, + h p + K ^= j 

x - cTj) ' v p (r) + i(l + xv) x • u p (r)) - g u [hi + hi + 
x ((a-i - a-j) • v w (r) + i(l + xs) (o-i x o-j) • u w (r)) - - 



+ r,- 



^ - Tj 



- ± Y A ) z (<r l + <7j)-Mr) (2) 

where M is the nucleon mass, v m (r) = {p,/ m (r)}, u m (r) = [p,/ m (r)], / m (r) = 
exp(— m m r)/47rr (with m m = m w ). The strong interaction parameters used by 
DDH are g nNN = 13.45, g pNN = 2.79, g uNN = 8.37, xv = 3.7, and xs = -0.12. The 
formula (2) is still used, [10,17]- [19] with these original strong interaction parameters. 

The effect we wish to incorporate is that Eq. (2) is not the complete PNC interac- 
tion between two-nucleons. This is because the PNC potential acting once, occurs in 
the midst of all orders of the strong potential, V. The resulting PNC reaction matrix 
G PNC (E) is a generalization of the Bruckner reaction matrix and is given by: 

G PNC (£) = &(E)V™ c n(E), (3) 
m) = 1 + -^--VQ(E) = 1 + T Ar r G(E), (4) 

hi 11Q III 11Q 

where E is an energy to be discussed below, and Q is an operator which projects onto 
unoccupied states. 

Our focus is on a first estimate of the effects of spin-isospin correlations. Thus we 
consider a nucleon of momentum P outside a core which is approximated as infinite 
nuclear matter. In this case, the PV single particle potential f/ PNC is given by: 

(M|f) PNC (P)| S ',t'> = ]T (P,s,t;k,m,t\G PNC (E)\P,s',t';k,m,m t ) A . (5) 

k,(fc<fcj?)m,mt 

The goal of this work is to assess the influence of tensor and spin dependent cor- 
relation effects on calculations of PNC observables. Since these are usually ignored, 
performing a schematic calculation seems appropriate. It might seem easiest to pa- 
rameterize the different contributions to Q(E) as function of r, the procedure of Miller 
& Spencer. Such a strategy will not work here, with the emphasis on finding hitherto 
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neglected contributions to the direct term. To see this consider the terms of the DDH 
potential of the form [p, f m ] ~ cr ■ r. The Miller-Spencer procedure would be to treat 
(r\Q(E)\K,,m,m t ) as a function of r, so that the direct matrix element of Eq. (5) 
would involve a vanishing volume integral of cr • r times a function of r. Indeed, the 
PV effective potential arises from the dependence of Q(E) on the relative momentum 

tz — (P — k)/2. (6) 

The importance of the spin-isospin correlations arises from the exchange of pions 
[11,12], so we separate the potential V into a one pion exchange term, Vope and a 
remainder approximated as being a central potential, V c . We include the effects of 
the central potential (which include the short-distance repulsion) to all orders and 
keeps terms of first-order in Vope- The application of the two-potential theorem to 
Eq. (4) gives 

n(E) » q c (e) + ni(E)—Q—v OP E n c (E), (7) 

hj — n 

with Q C (E) = 1 + E ® H V c Q C (E), and Q is the usual operator which projects onto 
two-particle-two-hole states. It is natural to model fl c (E) as 

(t\SI c (E)\k) = (1 + f c (r))e iK - r , (8) 

so that Eq. (7) can be expressed as 

(r\Q(E)\ K ) = (l + /c (r))e^ r + ^(r), (9) 

with 

ftc(r) = (1 + / c (r)) / ( ^ 3p2 ^ 2 e- pr ( P | Wl + /c)|«)- (10) 

The notation |) denotes a spatial overlap, so that V>At(r) are operators in spin-isospin 
space. The parameter uj (u 2 /M = —E), will be chosen to reproduce the results of 
Ref. [11]. Using a negative value of E causes the two-particle-two-hole fluctuations 
to have a short range of order ~ 1/uj. This allows us to neglect the effects of the 
operator Q because terms involving 1 — Q can be regarded as correction terms, of 
higher order in the density than the terms we examine here. 

The main result of this analysis is the simple model (9). It is necessary to compute 
the matrix elements of the pair-correlation operators of Eq. (1) to see if the results of 
Ref. [11] can be reproduced with such a simple formula. Evaluating Eq. (1) keeping 
only two-nucleon correlations leads to the result: 

/d k 
jT^Oikf - «)/(«)(«, S, M s , T, M T \Mr)\K, S, M s , T, M T ) A , (11) 

s,m s ,i,m T \ Z7r ) 

Mr) = -^^(E) J2S(r-n 3 )T t , Q(E), (12) 

i<j 
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where /(«) = 1 - §(£) + |(^) 3 and \k, S, M S ,T, M T ) A = \k, S, M s ,T, M t ) - 

(— l) s+T \ — k, S, M S ,T, M T ). We take k F = 1.36 fm -1 and choose the previously 
determined [14] function f c (r) = —e~ ar (1 — (3r 2 ) with a = 1.1 fm -2 , (3 = 0.68 fm~ 2 . 
The value of uj that leads to reproducing the results of Ref. [11] is u = 2fm~ 1 , which 
is of the expected order ~ k F - As shown in Fig. 1, and the present results reproduce 
the qualitative features of the pair-correlation functions of Ref. [11]. Furthermore, the 
integrals of our pp are in quantitative agreement with those of Ref. [11], so the present 
model should provide an adequate first assessment of the impact of such correlations 
on calculations of PNC observables. 
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FIG. 1. Pair-correlation functions. Solid curves are from Akmal & Pandharipande [11], 
dashed curves are from this work. 

Given our model of Eq. (9), we can now evaluate the operator f/ PNC of Eq. (5). 
Using Eq. (9) in Eq. (3), and keeping terms of first order in Pope leads to 

g pnc (£) » p pnc + ni(E)v OPE -®—v™ c + v™ c -Q—v OPE n c (E), (is) 

(14) 



= P PNC + AlT nu (£), 



PNC/ 
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where V™ c = (1 + f c )V PNC (l + f c ). The use of Eq. (13) in Eq. (5) specifies our model 
[16]. The numerical evaluations are straightforward, so we present the results. 

The operator £7 PNC (P) of Eq. (5) must be a pseudo-scalar spin-isospin operator. 
For our nuclear matter problem it is of the form 

U PNC (P) = ( T-P^( go + gi r z ). (15) 

For 133 Cs and 205 Th, (N - Z)/(N + Z) w 0.17, and we find that keeping the term 
gi amounts to keeping a correction to a term that is not large. So we take N = Z 
causing g x to vanish. To assess the importance of the spin-isospin correlation, we 
compare the influence of the two terms of Eq. (14) in Eq. (5). Denoting the results of 
using the first (\/ PNC ) term as and those of using the second term as Ago, with 
go = g^ + Ag , and evaluating the matrix element numerically leads to the results: 

g° = 24 U ~ 6-9 h° p - 3.6 h\ - 4.1 h° u - 4.1 hi (16) 
Ago = 2.43 f n - 1.2 h° p - .66 h\ - .351 h° u - .15 h x u (17) 

We see that the coefficient of each term is enhanced by about 10% if Ag is included. 
This is in contrast with many other nuclear structure effects which reduce the effects 
of PNC [8,9]. The present results constitute an argument that the effects of spin- 
isospin correlations need to be included in quantitative calculations of PNC effects, 
but are not large enough to not have a major impact [20]. Indeed, uncertainties due 
to nuclear structure effects might not have a large impact on PNC observables. 

But there is a more obvious source of change to present evaluations. Examination 
of the DDH potential, Eq. (2) immediately reveals the dependence on the product 
of strong and weak coupling constants. In a one-boson-exchange model one uses one 
strong and one weak meson-nucleon vertex function to construct the potential. But 
the strong coupling constants can be separately determined by computing the poten- 
tial V and choosing the parameters to reproduce experimentally measured scattering 
observables. This is the procedure of e. g. the Bonn potential [23]. As shown in 
Table I, the strong coupling constants required to fit phase shifts are much larger 
than those used originally by DDH (and used presently in Refs. [17-19]). 





DDH 


Bonn (OBEPR) 


9-kNN 


13.45 


13.68 


9p 


2.79 


3.46 


XV 


3.7 


6.1 


9u> 


8.37 


20 


XS 


-0.12 


0.0 



TABLE I. Comparison of strong coupling constants. 
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The coordinate-space potential (Table 14) of Ref. [23] is used for this comparison 
because it is a non-relativistic, local potential that is technically compatible with 
Eq. (2). The large values of xv an d g u are essential requirements to reproduce data, 
if a one-boson-exchange model is used [24]. 

Now consider the impact of using the Bonn parameters of Table I, in the DDH 
potential. If one considers a nucleon outside an N=Z core, the parameter combina- 
tions h° p ,l g p (l + xv) an d h^f-g^ determine the vector meson contributions to f/ PNC 
[5]. Thus using the Bonn parameters is equivalent to increasing the coefficient of the 
terms proportional to by a factor of 1.9. The effects of using values such as those 
of Table I have been noticed previously [25,26]. 
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FIG. 2. Constraints on the PNC coupling constants (xlO 7 ) that follow from using the 
Bonn coupling constants of Table I. 

To make a first assessment of the impact of this finding, we revise the results of Ta- 
ble VII and Fig. 9 of Ref. [17] by multiplying the coefficient of the terms proportional 
to h^ p by a factor of 1.9 [21]. This leads to Fig. 2. The allowed regions are between 
the dashed lines (pp), each set of solid lines ( 18 ' 19 F), and each set of dot-dashed lines 
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133 Cs, 205 T1. The results of the pp, and 18>19 F experiments can be explained with one 
set of weak-coupling parameters. The 205 T1 result with its large error band is almost 
consistent with that set of parameters, but the 133 Cs, and 205 T1 experimental results 
are not compatible, as pointed out earlier [17]. The results from compound nuclei are 
typically not shown in plots such as Fig. 2, but using the Bonn parameters of Table 
I, would reduce the discrepancies between theory and experiment. 

The constraints arising from (p, a) scattering experiments are not shown in Fig. 2, 
because currently interpreted provides constraints very similar to that of the 19 F 
experiment, and the calculations are not not complete. For example, U PNC , which 
should be complex, is treated as real. But it is worth remarking that an (n, a) 
experiment would provide very different constraints. The parity violation predicted 
using the original DDH potential is very small due to a cancellation between the 
pionic and vector meson exchange terms [22]. The use of the Bonn coupling constants 
(Table I) leads to a value of g n = 1.9 (see Eq. (18) of [13]) instead of g n = 0.2. This 
would be increased further by using a smaller value of f w (indicated in Fig. 2) than 
the DDH "best guess value" of Ref. [13]. If / w = 1, then g n = 3.2. 

The ranges of weak-coupling constants covered by Fig. 2 are roughly consistent 
with the DDH "reasonable ranges" and the same coupling constants account for much 
of the data, except for that of 133 Cs. This constitutes some success in explaining PNC 
phenomena, but the main improvement obtained here by using larger strong coupling 
constants could itself be true only within one-boson exchange models. All such models 
have large values of xv an d g^NN, but two-boson exchange models have smaller values 
of gu>NN [23]. Furthermore, the one-boson-exchange approach is not consistent with 
the present state of the art treatments of nucleon-nucleon scattering, see Refs. [27,28], 
and [29] . An appropriate treatment of PNC effects is one which involves incorporating 
PNC effects within an updated treatment of the nucleon-nucleon interaction. For 
example, Holstein has recently advocated an effective field theory treatment [30]. 

Thus our summary is that the present results demonstrate the need for an im- 
proved, updated, consistent incorporation of strong and weak interaction effects, and 
also indicate that the nuclear structure uncertainties might not be very severe. Thus 
the improvement of calculations of nuclear PNC effects seems to be an interesting 
and feasible task, even though much remains to be done. 
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